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Abstract 

two  part  study  was  undertaken.  In  the  first  part,  wlnglet 
effects  were  examined  on  a variety  of  wing  pi an forms.  A nonplanar 

lifting  surface  computer  program  was  used  to  calculate  aerodynamic 

• * 

coefficients  Including  lift,  induced  drag,  wing  pitching  moment, 
and  wing  root  bending  moment  coefficients.  Typical  cruise  flight 
conditions  were  examined  and  the  coefficients  calculated  for  the 
wing  alone  were  compared  to  those  obtained  for  the  wing  with  winglets 
Installed.  The  calculations  were  then  repeated  as  the  wing  aspect 
ratio,  sweep  angle,  and  dihedral  angle  were  varied.  Wlnglet  size  and 
orientation  remained  constant.  The  percentage  Induced  drag  reduction 
was  found  to  be  the  greatest  on  the  wing  with  the  highest  aspect  ratio  and 
wing  sweep  angle  (27. 5Z  for  an  aspect  ratio  7,  45  degree  swept  wing). 
Percentage  Induced  drag  reduction  Increased  slightly  with  Increasing 
positive  dihedral  angle.  The  greatest  Incremental  drag  reduction 
occurred  with  the  lowest  aspect  ratio,  highest  swept  wing. 

The  second  study  examined  the  effects  of  wlnglet  cant  angle  on 
induced  drag  reduction.  Vinglet  cant  angle  on  two  separate  wings 
was  varied  from  -3  degrees  to  +3  degrees  in  an  attempt  to  find  an 
optimum  value.  The  model  used  in  this  analysis  ignored  viscous 
effects  and  indicated  that  positive  cant  angles  (leading  edge  Inboard) 
yielded  the  greatest  reductions  in  induced  drag. 
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A PARAMETRIC  ANALYSIS  OF 
WINGLET  EFFECTS 


I.  Introduction 


I 


Objective 

The  objective  of  this  thesis  is  to  evaluate  the  effects  and 
potential  benefits  of  vlnglets  Installed  on  a variety  of  wing  planforms. 

A parametric  analysis  of  wlnglet  effects  will  be  conducted  to  determine 
the  possible  net  induced  drag  reduction  for  several  representative 
planforms.  The  planform  parameters  to  be  varied  include  wing  sweep 
angle,  aspect  ratio,  and  wing  dihedral.  The  analysis  will  be  conducted 
for  typical  cruise  conditions  only.  An  analysis  of  the  effect  of  wlnglet 
incidence  angle  will  also  be  made. 


Background 

■m 

The  present  concern  over  decreasing  energy  supplies  has  necessitated 
a critical  evaluation  of  Increasing  fuel  consumption.  Over  the  past 

j 

twenty-five  years.  United  States  energy  consumption  has  more  than 
doubled  (Ref  1:19).  The  National  Aeronautics  and  Space  Administration 
(NASA)  was  recently  tasked  by  Congress  to  establish  a comprehensive 
program  for  developing  aeronautical  fuel-conservation  technology.  In 
response,  NASA  established  a task  force  of  government  scientists  and 
NASA  engineers  to  obtain  recommendations  from  NASA  research  centers, 
aircraft  manufacturers,  and  other  government  agencies.  This  task 
force  has  considered  numerous  improvements  in  all  areas  of  aeronautical 
technology.  One  aerodynamic  improvement  under  investigation  is  the  use 
of  winglets  (Ref  1:20). 
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The  United  States  Air  Force  (USAF)  is  concerned  with  increasing 
fuel  costs  and  the  threat  of  future  Increases.  The  USAF  is  evaluating 
the  use  of  winglets  on  current  USAF  aircraft.  The  KC-135  and  the 
C-141,  which  are  two  heavy  fuel  consumers » are  prime  candidates  for 
modification  with  winglets. 

Winglets 

Winglets,  as  used  in  this  study,  are  small,  nearly  vertical  wing- 
like  surfaces  mounted  at  the  tips  of  a wing.  These  devices  are  mounted 
nearly  perpendicular  to  the  main  wing  and  are  normally  canted  slightly 
outward.  They  act  as  vortex  diffusers  and  provide  reductions  in  induced 
drag  greater  than  those  achieved  with  simple  wing  tip  extensions 
(Ref  2:1). 

For  many  years  it  has  been  recognized  that  a nonplanar  lifting 
system  should  have  less  Induced  drag  than  a planar  wing  (Ref  3:1).  As 
early  as  1897  a patent  was  obtained  by  Lanchester  for  vertical  surfaces 
at  the  wing  tips  (Ref  3:1).  Since  that  time  a number  of  theoretical 
analyses  have  indicated  the  significant  improvements  possible  with  non- 
planar systems  (Ref  3:2). 

In  addition  to  the  reduction  in  Induced  drag,  winglets  produce 
several  other  effects  which  relate  to  the  total  drag  improvement.  First, 
the  profile  drag  of  the  winglet  itself  negates  some  of  the  drag  improve- 
ment. Second,  for  most  wing  and  winglet  combinations,  the  required  lift 
will  be  produced  at  a lower  angle  of  attack.  Thus,  the  parasite  drag  of 
the  wing  will  be  decreased.  Third,  at  the  intersection  of  the  wing  and 
winglet,  interference  effects  cause  a slight  increase  in  drag. 


1 

i 


- 

i 

i 

\ 

| 


i « 


HE 


, * 


l : 
■ 


2 


Ocher  effects  which  must  be  considered  are  the  possible  weight 
increases  and  trim  penalties.  In  addition  to  the  weight  of  the  wlngleta 
and  attachments,  increased  aerodynamic  loads  on  the  main  wing  may  require 
strengthening  the  wing  and  increasing  its  weight.  Part  of  the  vinglet 
benefits  may  be  negated  due  to  trim  effects.  As  the  loading  at  the 
wing  tip  is  increased,  a net  nose  down  pitching  moment  may  occur.  This 
must  be  balanced  by  a moment  in  the  opposite  direction  and  may  require 
a down-load,  or  lift  reduction  on  the  tall  of  the  aircraft.  Thus, 
potential  benefits  of  wlnglets  must  be  weighed  against  their  penalties, 

I 

and  care  must  be  taken  to  look  at  net  improvement  only. 

Previous  Investigations 

Dr.  Richard  R.  Whitcomb  of  the  NASA  Langley  Research  Facility  has 
conducted  experiments  to  determine  the  effects  of  wlnglets  on  overall 
aircraft  performance  (Ref  3) . His  studies  have  dealt  primarily  with 
the  potential  benefits  from  winglet  addition  to  first  and  second  genera- 
tion jet  transport  wings. 

The  Boeing  Commercial  Aircraft  Company  has  recently  completed  a 
study  of  the  effects  produced  by  wlnglets  on  the  KC-135  and  C-141 
(Ref  4).  A winglet  parameter  study  was  also  conducted  by  Boeing  to 
determine  the  potential  aerodynamic  improvement  of  the  KC-135  (Ref  4: 
108-169).  The  parameters  varied  were  leading  edge  sweep  angle,  taper 
ratio,  area,  length,  and  winglet  dihedral  angle.  In  all  cases  the 
wlnglets  were  located  on  the  upper  surface  of  the  main  wing.  An  aero- 
dynamic vortex  lattice  computer  program  was  used  in  the  Boeing  study. 

Investigations  to  date  have  dealt  primarily  with  the  addition  of 

( 

wlnglets  to  transport  aircraft  such  as  the  C-141  and  the  KC-135.  The 
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parametric  analysis  of  winglets  with  planform  parameters  described 
previously  has  yet  to  be  reported.  This  parametric  Investigation  was, 
therefore,  considered  to  be  a fruitful  area  of  research. 


Scope 

This  analysis  Is  limited  to  the  investigation  of  wlnglet  effects 
on  trapezoidal  planforms  at  typical  cruise  flight  conditions.  Only  the 
effects  of  upper  surface  winglets  will  be  examined.  The  size,  location, 
and  orientation  of  these  winglets  will  remain  constant,  except  during 
the  wlnglet  cant  angle  study. 


Approach 

To  conduct  this  analysis,  a new  nonplanar  lifting  surface  theory  is 
used.  This  method,  developed  by  Mr.  M.  I.  Goldhammer  of  Douglas  Aircraft 
Company,  represents  nonplanar  lifting  systems  by  distributed  vorticity 
(Ref  5 and  6) . This  method  has  been  used  to  compute  the  induced  drag 
benefit  of  winglets  and  the  agreement  of  computed  results  with  NASA 
experimental  data  is  excellent  (Ref  5:51-61).  The  validity  of  this 
method  will  be  further  examined  in  the  next  chapter. 

The  Nonplanar  Lifting  Systems  (NPLS)  method  is  a nonlinear,  non- 
planar, lifting  surface  theory.  The  term  lifting  surface  theory 
indicates  that  wing  thickness  effects  are  neglected.  Wing-like  bodies 
are  represented  by  a cambered  surface  with  zero  thickness.  This  theory 
makes  no  small  perturbation  assumptions  as  do  previous  lifting  surface 
theories  (Ref  5:5). 

The  basic  problem  is  to  determine  aerodynamic  forces  and  moments 
on  an  arbitrary  three  dimensional  lifting  system.  The  system  is  assumed 
to  he  placed  In  a uniform,  inviscid,  and  irrotntional  flow  and  is 
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represented  by  an  array  of  rectangular  elements.  The  calculation  of 


aerodynamic  forces  and  moments  Is  accomplished  by  the  digital  computer. 

A computer  program  utilising  the  NPLS  theory  has  been  written  by  Douglas 
Aircraft  Company  and  Is  employed  in  this  thesis.  This  program  gives 
aerodynamic  force  and  moment  coefficients  for  either  a single  wing  or 
a combination  of  lifting  surfaces.  By  calculating  the  aerodynamic 
forces  and  moments  for  the  wing-wlnglet  combination  and  comparing  them 
to  those  of  the  wing  alone,  the  Invlsdd  Incremental  effects  of  the 
wlnglets  are  determined.  A more  detailed  summary  of  the  NPLS  method 
Is  contained  In  Appendix  A. 
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II.  Use  and  Verification  of  the  NPLS  Method 

Nonplanar  Lifting  Systems  Program 

The  basic  tool  used  to  conduct  this  investigation  was  the  Douglas 
Aircraft  Company  Nonplanar  Lifting  Systems  (NPLS)  Program.  Detailed 
instructions  for  the  use  of  this  program  may  be  found  in  Reference  5. 

The  analytical  method  used  in  this  program  was  a finite  element 
technique  which  required  the  numerical  definition  of  the  lifting  system 
to  be  evaluated.  Since  a lifting  system  was  represented  by  an  array  of 
rectangular  elements,  the  geometric  characteristics  of  this  array  .were 
first  calculated  and  then  input  to  the  computer  program. 


Program  Input 

To  analyze  a particular  configuration  using  this  program  it  was 
necessary  to  geometrically  define  each  surface  or  combination  of  sur- 
faces. Component  surfaces  of  a system  were  separated  and  each  surface 
was  treated  Independently.  Each  surface  was  initially  orientated  in  a 
defining  (x,  y,  z)  system  and  then  rotated  into  the  desired  position. 

To  define  a given  planform,  the  program  required  that  the  leading 
and  trailing  edge  coordinates  at  the  tip,  root,  and  either  side  of  a 
discontinuity  be  given.  These  required  points  are  illustrated  in  Fig.  1. 

Once  the  planforms  were  defined,  each  surface  was  divided  into 
spanwise  sections.  This  division  was  based  on  two  factors.  First,  since 
the  planform  geometry  was  represented  by  rectangular  elements,  the  span- 
wise  divisions  needed  to  adequately  represent  the  resulting  "stairstep" 
leading  and  trailing  edges.  The  second  factor  to  be  considered  was  the 
anticipated  gradients  in  spanwise  loading.  Since  this  method  assumed 
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that  spanwlse  loading  was  constant  over  each  section,  close  spacing  was 
used  In  regions  where  large  load  variations  were  anticipated.  It  was 
suggested  that  section  widths  be  limited  to  2-5Z  of  the  wing  semi-span 
at  the  wing  tip  (Ref  5:67).  To  adequately  represent  a 45  degree  swept 
wing,  14  divisions  were  used.  This  representation  with  the  location 
of  these  spanwlse  divisions  Is  shown  in  Fig.  2. 

Once  the  spanwlse  divisions  were  selected,  the  chordwlse  divisions 
were  chosen  such  that  the  complete  planform  was  represented  by  an  array 
of  rectangular  elements.  Chordwlse  divisions  were  also  made  such  that 
closely  spaced  elements  were  used  in  regions  where  rapid  changes  in 
loading  were  expected. 

The  Elementary  Vortex  Distribution  approach  used  In  this  method 
assumes  that  the  functional  form  of  the  vorticlty  distribution  on  the 
leading  edge  element  is  F(l/*^j),  where  s is  a streamwise  coordinate. 

On  the  remaining  elements  a trapezoidal  distribution,  F(s),  is  assumed 
(Ref  5:18-23).  This  distribution  results  from  the  superposition  of  two 
triangular  distributions  (Fig.  A-l).  Since  the  loading  distribution  on 
the  leading  edge  is  well  represented  by  the  inverse  square  root  vortex 
distribution.  It  is  not  necessary  to  make  the  leading  edge  element  very 
small.  For  this  reason,  the  leading  edge  element  was  allowed  to  be 
rather  large  (5  per  cent  of  chord)  in  relation  to  the  remaining  elements. 
The  second  element  was  chosen  to  be  smaller  (2-1/2  per  cent  of  chord) , 
which  allowed  a smooth  transition  to  the  triangular  distribution.  Twelve 
chordwlse  divisions  were  used  in  this  investigation.  These  divisions 
and  their  chordwlse  locations,  are  illustrated  on  the  root  section  in 
Fit*.  2.  A representation  of  this  type  was  recommended  by  Douglas  (Ref  5: 


The  dihedral  angle  of  a wing,  T,  can  be  specified  in  either  of  two 


! 
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way 8.  It  can  be  specified  by  inputting  the  vertical  displacements  of 
the  leading  edge  from  the  x-y  plane  or  the  wing  can  be  initially  speci- 
fied to  be  entirely  in  the  x-y  plane  and  then  rotated  through  the 
required  dihedral  angle.  The  first  method  is  the  most  accurate  since 
a rotation  through  an  angle,  T,  reduces  the  span  by  a factor  of  (l-cosO 
and  decreases  the  reference  area.  This  would  cause  a 1.5  per  cent 
decrease  in  reference  area  for  a dihedral  angle  of  10  degrees. 

Winglet  Definition 

The  winglets  were  defined  using  the  above  procedure.  Seven  span- 
wise  and  twelve  chordwlse  divisions  were  used  to  fully  define  each 
winglet  in  this  study.  The  winglet  is  initially  defined  in  an  (x,  y,  z) 
system.  It  is  then  translated  and  rotated  to  join  with  the  wing  at  the 
wing  tip.  The  cant  angle  rotation  is  done  about  the  y axis  of  the 
defining  system;  the  dihedral  angle  rotation  is  done  about  the  x axis. 
The  surface  is  then  translated  to  its  final  position  by  specifying  the 
coordinates  of  the  winglet  apex  (x^,  yQ,  z^) . This  procedure  is  illus- 
trated in  Fig.  3. 

Wlng-Winglet  Junction 

The  junction  of  wing  and  winglet  had  to  be  completed  such  that  gaps 
were  avoided  and  any  overhang  was  minimized  (Ref  5:69-72).  The  widths 
of  the  wing  rectangular  elements  at  the  wing  tip  were  selected  to  be 
large  enough,  (4  per  cent)  to  accommodate  the  expected  ranges  of  inci- 
dence, + 5 degrees,  and  camber  of  the  winglet.  The  width  of  the 
rectangular  elements  at  the  winglet  root  was  15  per  cent  of  the  winglet 
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span  to  allow  for  the  wing  camber.  This  prevented  Intersecting  or 
passing  near  boundary  condition  points  on  the  other  surface. 


Program  Output 

The  output  of  this  program  Includes  the  basic  Input  parameters , 
detailed  geometric  data,  distributions  of  pressure  jump  coefficients 
on  each  element,  and  aerodynamic  loading  data.  This  data  Is  given  for 
each  system  surface  and  for  the  complete  system  at  each  angle  of  attack. 

The  aerodynamic  coefficients  calculated  by  this  program  are  the 
lift  coefficient  (C^),  the  Induced  drag  coefficient  (C^),  the  side  force 
coefficient  (C^) , and  the  pitching  (C^ , rolling  (CR) , and  yawing  (C^) 
moment  coefficients.  For  symmetric  configurations,  total  values  of  C^, 
C„,  and  C_  are  for  the  right  half  of  the  system  only.  The  actual  total 
values  are  zero.  The  total  pitching  moment  is  referenced  about  the 


quarter  chord  of  the  mean  aerodynamic  chord.  Yawing  and  rolling  moments 
are  referenced  about  the  origin  (Fig.  3).  These  surface  total  aerodyna- 
mic coefficients  are  normalized  by  surface  area,  span,  and  reference 


chord. 


Convergence 


The  convergence  characteristics  of  this  solution  have  been  shown 
to  be  quite  good  (Ref  5:32-36).  Values  of  C are  equally  sensitive  to 

L 

the  number  of  chordwise  and  spanwise  divisions.  Values  of  are  most 

sensitive  to  the  number  of  chordwise  divisions.  Values  of  C are  most 

II 

sensitive  to  the  chordwise  arrangement  of  elements  (Ref  5:33). 

This  program  will  accept  as  many  as  500  finite  elements;  252  were 
used  to  define  a wing-winglet  combination.  Studies  completed  by  Douglas 
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indicate  that  10  spanwise  and  chordwise  divisions  yield  lift  coefficients 
within  1.5  per  cent,  drag  coefficients  within  4.0  per  cent,  and  moment 
coefficients  within  5 per  cent  of  those  obtained  with  a 400  element 
solution  (Ref  5:34).  Computational  time  increases  greatly  with  the 
number  of  elements  used.  It  is  felt  that  the  number  selected  for  this 
investigation  represents  each  planform  accurately  and  achieves  a 
desirable  balance  between  accuracy  and  computation  time. 


Method  Verification 

The  Douglas  Aircraft  Company  has  compared  results  from  the  NPLS 
method  both  with  experimental  results  and  with  other  theoretical  methods. 
(Ref  5:32-62).  Calculation  of  induced  drag  reduction,  using  this  method 
agrees  with  NASA  experimental  results  for  the  KC-135A  within  11  per  cent 
(Ref  5:59).  Flow  interference,  viscous,  and  aeroelastic  effects  not 
considered  in  this  method  may  account  for  this  difference. 

One  test  Douglas  conducted  of  the  NPLS  method  was  to  calculate  the 
aerodynamic  coefficients  of  the  aspect  ratio  4,  rectangular  wing  dis- 
cussed in  Ref  7.  The  author  attempted  to  duplicate  the  Douglas  results 
to  verify  understanding  and  successful  application  of  the  NPLS  method. 
Experimental  data  for  this  wing  with  a variety  of  endplates  (winglets) 
was  available  as  was  data  from  Douglas  NPLS  computations.  This  configura- 
tion used  a NACA  64^A412  airfoil  section.  The  ordinates  of  this  section 
are  listed  in  Table  I. 

The  program  input  was  completed  using  the  technique  described 
previously.  Although  fewer  rectangular  elements  could  have  easily 
represented  this  rectangular  planform,  168  elements  were  used  (14  span- 
wise  X 12  chordwise).  This  representation  was  selected  since  it  was  to 
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be  used  for  the  wings  in  the  remaining  representations. 

The  planform  parameters  were  non-dimensionalized  by  the  wing 
semispan  (b/2)  and  are  as  follows 


b - 2 


AR  - 4 


S - 1 


c ■ .5 


The  leading  and  trailing  edge  coordinates  necessary  to  define  this 
planform  are  shown  in  Fig.  4.  The  camberline  was  approximated  by  the 
airfoil  meanline,  the  average  of  the  upper  and  lower  airfoil  ordinates 


"c  * 2 ("l  + V 


where  n represents  the  y ordinate  in  percent  chord.  The  subscripts  c, 
1,  and  u refer  to  the  camberline,  upper  airfoil  surface  and  lower 
airfoil  surface. 

The  camber  angle  at  the  midpoint  of  each  rectangular  element  was 


calculated  analytically  from 


-1 


e “ tan 


In  Eq  (3),  At^  represents  the  change  in  camberline  and  Ax  repre- 
sents the  change  in  the  x ordinate  of  the  rectangular  element,  c 
represents  the  camber  angle  at  the  midpoint  of  the  rectangular  element. 
This  camber  angle  is  positive  in  the  angle  of  attack  sense.  Calculated 
values  for  this  airfoil  are  listed  in  Table  II.  A comparison  of  the 
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NPLS  approximation  and  the  actual  meanline  is  shown  in  Fig.  5. 

The  flight  condition  investigated  was  Mach  ■ .211.  Angles  of 
attack  ranged  from  zero  through  ten  degrees. 

The  results  calculated  by  the  author  are  shown  in  Fig.  6 and  are 
compared  to  the  experimental  results.  As  can  be  seen  from  this  compari- 
son, both  CLQ  and  are  very  well  predicted.  The  shape  of  the  drag 
polars  agrees  very  well  up  to  - .8  where  experimental  results 
indicate  flow  separation  begins.  The  magnitude  differences  in  these 
polars  is  due  to  the  NPLS  prediction  of  only  induced  drag.  The  slope 
of  the  pitching  moment  curve  up  to  flow  separation  is  well  predicted 

although  C is  slightly  underpredicted.  This  discrepancy  is  consistent 
MO 

with  other  lifting  surface  theories  and  is  probably  due  to  thickness 
effects  (Ref  5).  These  results  are  consistent  with  Douglas's  calcula- 
tions. 

Figure  7 shows  a comparison  with  experimental  data  for  the  same 
airfoil  with  a rectangular  winglet  attached.  This  winglet  was  a full 
chord  winglet  attached  at  a 90  degree  dihedral  angle  and  a zero  degree 
cant  angle.  The  winglet  height  was  .765  of  the  main  wing  chord 
(.3825  b/2).  The  camber  distribution  used  for  the  winglet  was  that  of 
a flat  endplate. 

As  shown  in  Fig.  7,  and  C^a  are  very  well  predicted  and  the 

drag  prediction  is  in  very. good  agreement  through  e .8.  This  method 

fails  to  accurately  predict  aerodynamic  data  after  flow  separation 

begins.  is  now  slightly  over  predicted  while  the  slope  is  still 

MO 

in  very  close  agreement  with  experimental  data. 
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Aerodynamic  Characteristics  oi 
Wing  With  Endplate  (Calculated 


Fig.  6.  Aerodynamic  Characteristics  of  AR=4  Rectangular 
Wing  (Calculated  vs  Experimental) 


rodynamic  Characteristics  of  AR=4  Rectangular 
ng  With  Endplate  (Calculated  vs  Experimental 
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The  exact  representation  used  by  Douglas  to  aodel  this  wing  was 
not  presented  in  References  5 or  6.  The  author's  aodel  was  constructed 
Independently. 

A comparison  of  the  author's  calculations  with  those  calculated  by 
Douglas  (Ref  6:8)  shows  very  little  variation  In  lift  coefficients  and 
moment  coefficients  (2Z).  The  drag  coefficients  presented  In  Ref  6 
include  a skin  friction  model  and  were  a constant  .010  higher  than  the 
author's  for  the  wing  alone  and  .012  higher  for  the  wing  with  endplates 
attached. 

This  comparison  of  NPLS  calculations  with  experimental  data  served 
as  a validation  for  the  model  and  the  ability  of  the  author  to  use  It 
correctly.  The  next  chapter  describes  hew  the  model  was  used  In  the 
current  parametric  study. 
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III.  Calculation  of  Aerodynamic  Coefficients 
Wing  Parameter  Study 

The  object  of  this  Investigation  was  to  determine  the  effects  of 
vinglet8  on  a variety  of  planforms.  Wlnglets  have  been  Investigated 
primarily  with  regard  to  their  installation  on  transport  wings  similar 
to  that  of  the  KC-135. 

This  thesis  examines  winglet  effectiveness  on  a variety  of  plan- 
torms  in  order  to  evaluate  the  benefits  which  could  be  expected  on 


other  types  of  aircraft  and  provide  guidance  for  further  study.  The 
wing  parameters  which  were  varied  were  aspect  ratio,  wing  sweep  angle, 
and  dihedral  angle.  The  variation  of  these  parameters  and  comparison 
of  the  resulting  aerodynamic  data  offers  insight  into  winglet  effects 
on  a variety  of  military  and  civil  aircraft. 

Although  the  NPLS  method  allows  the  representation  of  rather  complex 
planforms,  simple  trapezoidal  planforms  were  selected  for  this  study.  A 
tabular  listing  of  the  wing  parameters  which  were  varied  is  given  in 
Table  III.  These  wings  are  also  illustrated  in  Fig.  8. 


Winglet  Cant  Angle  Study 

The  theory  and  experimental  results  of  Reference  3 indicate  that  in 
order  to  be  effective,  winglets  must  produce  significant  side  forces  to 
reduce  the  lift  induced  inflow  at  the  wing  tip.  NASA  results  indicate 
that  the  optimum  winglet  effectiveness  occurs  at  negative  (leading  edge 
outboard)  cant  angles  (Ref  3).  A study  by  Douglas  Aircraft  Company 
indicates  that  at  positive  cant  angles  the  greatest  induced  drag  reduc- 
tion occurs  (Ref  5:45). 
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To  further  Investigate  the  effects  of  vinglet  cant  angle  on  induced 
drag  reduction,  the  cant  angle  on  two  planforos  was  varied  between  -3 
degrees  and  +3  degrees  in  an  attempt  to  find  an  optimum  value. 

Program  Input 

The  procedure  used  in  both  of  the  above  Investigations  was  identi- 
cal. In  each  case  it  was  necessary  to  represent  the  lifting  surface,  or 
lifting  system,  as  an  array  of  rectangular  elements.  Once  the  necessary 
divisions  were  made,  and  the  representation  was  complete,  it  remained  to 
input  the  necessary  wing  and  winglet  parameters  and  definition  points  to 
the  NPLS  Computer  Program  to  calculate  the  required  aerodynamic  coeffl- 

I 

cients.  Aerodynamic  data  was  calculated  for  five  angles-of-attack  with 
winglets  and  compared  to  that  for  the  wing  alone  to  determine  the 
incremental  effects  due  to  winglet  addition. 

Airfoil 

The  basic  airfoil  section  used  to  construct  the  wings  in  this 
investigation  was  the  Boeing  BAC  311  airfoil  section.  The  camberline 
was  approximated  by  the  meanline  of  the  airfoil  ordinates  with  the 
following  relation 

"c  ■ i <"i + V <2> 

Twelve  chordwise  elements  were  used  to  obtain  the  approximation 
shown  in  Fig.  9.  Table  IV  lists  the  slopes  of  the  camberline  (in  degrees) 
which  were  input  to  the  NPLS  program.  These  values  were  determined  analy- 
tically from 

-1  Anc 

e = tan  -r-  (3) 

Ax 


* f i 
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In  Eq  (3),  e represents  the  camber  angle  at  the  midpoint  of  each  rectan- 
gular element,  x is  the  chordwlse  location  measured  in  precent  of  total 
chord. 


Planform  Definition 

Fourteen  spanwlse  divisions  were  made  using  the  criteria  described 
in  Chapter  II.  Thus,  each  wing  was  represented  by  168  (12  X 14) 
rectangular  elements. 

For  simplicity,  all  wing  distance  parameters  were  non-dimensionalized 
by  the  wing  semispan  (b/2).  Since  all  configurations  analyzed  were  sym- 
metric, only  the  right  half  of  each  planform  was  defined.  The  NPLS 
program  formed  the  entire  wing  using  a mirror  image  technique.  Given 
the  aspect  ratio,  AR,  the  non-dimensional  reference  area,  S,  was  found 
from  the  non-dimensional  wing  span,  b,  using  the  definition  of  aspect 
ratio 

2 

« ' T w 


The  taper  ratio,  X,  remained  constant  (X  = 1/3)  for  each  planform.  The 
root  chord  and  tip  chord  length  were  found  from  the  following  relation 
for  trapezoidal  wings. 


AR 


2b 

cr  (1  + X) 


ct  = xcr 


^‘l/4c  Was  sPec*^*e^  f°r  each  planform.  The  leading  and 
coordinates  required  to  define  each  planform  were  found 


(5) 


trailing  edge 
using 


tan  A 


LE 


= tan  A 


1/Ac 


J.  0 - X) 
AR  (1  + X) 


(6) 
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1 Meanline 


tanA  . 3_  1L=_*1 

CanAl/4c  AR  (1  + X) 


(6) 


tanA 


TE 


The  mean  aerodynamic  chord  and  location  of  the  moment  center  for  a 
trapezoidal  wing  are  given  by  the  following  relations  (Ref  8:3,  17) 


c 


2 ,1  + 1 + X2, 

3 ' 1 + X } 


m b (1  + 2X) 
ymc  " 6 (1  + X ) 


c 

x - + y tanA. 

me  4 me  l/4c 


(7) 


All  of  the  required  points  and  parameters  necessary  to  fully  define  a 

I 

wing  for  the  NPLS  program  are  Illustrated  in  Fig.  10.  The  parameters 
calculated  for  the  planforms  in  this  study  are  listed  in  Table  V. 

Winglet s 

The  winglet  selected  for  use  in  this  investigation  was  similar  to 
that  selected  for  use  on  the  KC-135  (Ref  4:12).  Although  some  investi- 
gations (Ref  2 and  3)  to  date  have  examined  the  effects  of  both  an  upper 
and  a lower  winglet,  only  the  upper  winglet  was  examined  here.  This 
winglet  (Fig.  11)  was  mounted  on  all  planforms  with  an  incidence  angle 
of  zero  degrees  for  the  planform  analysis.  Although  NASA  tests  indicate 
a substantial  negative  incidence  angle  may  be  used,  zero  degrees  was 
selected  as  a baseline  value  for  this  study  (Ref  2 and  3).  In  all  cases 
the  winglet  was  mounted  with  the  winglet  trailing  edge  flush  with  the 
wing  trailing  edge  and  mounted  at  a 70  degree  dihedral  angle.  The  wing- 
let was  cambered  inboard  to  produce  the  Inward  loads  necessary  to  reduce 
the  lift  induced  inflow  and  decrease  the  induced  drag.  This  location 


at 

H 


« 


at 

H 

at 

A 

4-» 

3 

a 


g 

at 

ex 

o 

u 

a, 

CO 

>-l 

a, 

z 


w 

H 

< 


<D 

0 

W 2 

- 0 

^ >o 


o 

-*N* 

■< 


o 


CO 


a 

< 


ext 

c 


4te 


ojinor-oovOOLfi 

OHO'rtOCOHO 

■CO  to  0 H H SO  CJ 
H N M I I CO 


vOO'iouiia<HNx> 

•«fvomt^inoiosD 

r'NvOO'ON^CI 

o a t>‘  o O'  o h u 
~«rgco^*— ct'loi'l' 


minininxiinintn 


O' 

O' 

O' 

o 

r- 

o 

o 

O 

fM 

04 

04 

sO 

o 

o 

o 

•<* 

■<* 

■sO 

o 

in 

m 

— i 

fH 

H 

H 

04 

04 

04 

nOsOvOvOOOOO 

OOOOOOOOOOOO 

NNNNOOmin 


(OtOtOtOtOtOtOtO 

tototocototototo 

cotocototocotoco 


< — < *-H  I — I r—i  \Q  O O O 

r~c^r-c^voooo 

mioiniovooooo 


NMIMNNIMNN 


•o 

X 


XI 

x 


u 

g 

>> 


lo 


00 

■ S -V: 

t> 


'O’tNvOcjino'Co 

'H^ooTtNW'ts 

'tvOO'iMininom 


(OsCT'Tj'vOvOO 

no  ooiroorooom 

O'cOffl'tOO'hN 

iMmm'too'^vo 


rHOCT'OOOcOOO 

IflHHOOvOOWC' 

aiosointfvo®H 

(OinNonnnH 


O' 

in 

CO 

CO 

to 

'O 

vO 

O' 

r- 

n- 

o 

to 

oo 

'tf 

v£> 

o 

m 

F-< 

r~ 

r- 

<M 

O 

r—i 

H 

co 

o 

o 

o 

C4 

o 

o 

o 

o 

o 

o 

o 

O 

vOnD'DvDnOO'O'O 

H rH  c-H  r-H  H rH  H H 


oo-'t'O'oovOvO-Hr^ 

tOHOOnivOHOWl 

HOtM'JM'OO'Kt 

NtOOXlNNN^ 


inminuiHfor'N 

O'^O'O'HrtHH 

OOOOsOtOTf't 

cOCOrOfOtO'^miO 


28 


V y 


and  placement  has  been  shown  to  be  near  optimum  (Ref  3:6). 

The  vinglet  airfoil  section  used  in  this  analysis  is  the  NASA 
general  aviation  airfoil  described  in  Reference  3.  The  ordinates  of 
this  airfoil  are  listed  in  Table  VI.  The  camber  angles  for  this  air- 
foil were  calculated  as  described  previously  and  are  listed  in  Table  VII. 
Figure  12  shows  the  resulting  camberline  approximation  used  In  this 
analysis. 

Winglet  Definition 

Each  winglet  was  represented  by  84  (7  spanwise  and  12  chordwlse) 
rectangular  elements  based  on  the  criteria  covered  in  Chapter  II.  By 
inputting  the  winglet  aspect  ratio  (2.32),  taper  ratio  (.338)  quarter- 
chord  sweep  (32.918  degrees),  spanwise  divisions,  and  chordwlse 
divisions,  the  winglet  was  fully. defined  ./The  winglet  orientation  was 
specified  by  inputting  the  cant  and  dihedral  angles  illustrated  pre- 
viously in  Fig.  3.  Winglet  location  wa6  specified  by  inputting  the 
coordinates  (xQ,  yQ,  zQ)  of  the  winglet  apex.  The  yQ  and  zQ  coordinates 
were  identical  to  the  wing  tip  coordinates  when  winglet  cant  was  zero. 

The  value  of  xQ  was  calculated  which  placed  the  trailing  edge  of  the 
winglet  flush  with  the  trailing  edge  of  the  main  wing 


x_  * c , . . + y.  tan  A__  - c , . , (8) 

0 r (wing)  70  TE  r (winglet)  7 

When  a positive  winglet  cant  angle  was  used,  the  apex  of  the  winglet  was 
moved  inboard  to  prevent  a gap  between  the  wing  and  winglet.  The  new  y 
coordinate  of  the  winglet  apex  was 


Cr  (winglet)  s*na  cant 


(9) 
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Table  VI 


Airfoil  Coordinates  for  Winglets 


x/c 

y/c  for  1 

Upper 

Surface 

Lower 

Surface 

0 

0 

0 

.0020 

.0077 

-.0032  • 

.0050 

.0119 

-.0041 

.0125 

.0179 

-.0060 

.0250 

.0249 

-.0077 

.0375 

.0296 

-.0090 

.0500 

.0333 

-.0100 

.0750 

.0389 

-.0118 

. 1000 

.0433 

-.0132 

. 1250 

.0469 

-.0144 

. 1500 

.0499 

-.0154 

. 1750 

.0525 

-.0161  • 

.2000 

.0547 

-.0167 

.2500 

.0581 

-.0175 

.3000 

.0605 

-.0176 

.3500 

.0621 

-.0174 

.4000 

.0628 

-.0168 

.4500 

.0627 

-.0158 

.5000 

.0618 

-.0144 

.5500 

.0599 

-.0122 

. 5750 

.0587 

-.0106 

. 6000 

.0572 

-.0090 

.6250 

.0554 

-.0771 

.6500 

.0533 

-.0052 

.6750 

.0508 

-.0033 

.7000 

.0481 

-.0015 

.7250 

. 0451 

.0004 

.7500 

.0419 

.0020 

.7750 

.0384 

.0036 

.8000 

.0349 

.0049 

.8250 

.0311 

.0060 

.8500 

.0270 

.0065 

.8750 

.0228 

.0064 

. 9000 

.0184 

.0059 

.9250 

.0138 

.0045 

. 9500 

.0089 

.0021 

.9750 

.0038 

-.0013 

1.0000 

-. 0020 

-.  0067 
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Fig.  12.  NASA  Winglet  Airfoil  M< 


NPLS  Program  Use 


Once  the  required  parameters  were  calculated,  they  were  arranged 
in  the  proper  input  format  (Ref  5:72-80). 

The  flight  condition  was  specified  by  Inputting  the  desired  Mach 
number  and  as  many  as  five  separate  angles  of  attack.  Mach  .7  was 
selected  for  use  in  this  investigation.  Angles  of  attack  of  zero  through 
six  degrees  were  investigated.  These  safely  represent  those  flight  con- 
ditions for  which  flow  separation  is  not  a factor. 

The  computer  program  has  provisions  to  check  input  data  for 
inconsistencies,  however,  the  input  was  closely  rechecked  to  Insure  the 
desired  representation  was  obtained. 

For  the  wing-winglet  combinations  used  in  this  analysis  (252  rec- 
tangular elements) , approximately  130  CP  seconds  execution  time  were 
required. 


Program  Output 

The  program  output  for  each  case  began  with  a listing  of  the  basic 
geometric  parameters.  This  was  followed  by  a detailed  element-by-element 
listing  of  geometric  data.  This  data  was  carefully  compared  to  the  input 
data  to  insure  the  desired  representation  was  obtained.  Local  dihedral 
angles  were  carefully  checked  to  insure  the  surface  was  not  excessively 
warped.  Had  an  excessive  local  dihedral  angle  existed,  the  representa- 
tion would  have  been  altered  to  better  represent  the  desired  wing.  Mo 
automatic  check  for  surface  gaps  exists  in  this  program.  Only  be  care- 
fully analyzing  the  aerodynamic  data  could  an  irregularity  such  as  a gap 
between  wing  and  winglet  be  detected. 
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The  NPLS  program  does  not  include  a plotting  routine.  Since  all 
output  was  in  tabular  form,  a plotting  routine  was  written  to  display 
graphically  the  calculated  coefficients.  Samples  of  these  aerodynamic 
plots  are  found  in  the  next  chapter;  the  remainder  are  located  In 
Appendices  B and  C. 


f 

c 

J 


il 


IV.  Results 


Wing  Parameter  Study 

Figures  13-18  are  representative  of  the  aerodynamic  plots  contained 
In  Appendix  B.  These  plots  Illustrate  the  effect  of  vlnglet  addition  to 
each  particular  wing  Investigated.  These  curves  show  the  lift  coeffi- 
cient as  a function  of  angle  of  attack  (Figure  13) , Induced  dreg 
coefficient  (Fig  14),  pitching  moment  coefficient  (Fig.  16),  and  root 

bending  moment  coefficient  (Fig.  17).  Also  shown  is  lift  coefficient 
2 

squared,  C^»  as  a function  of  the  induced  drag  coefficient  (Fig.  15) 
and  the  variation  In  spanwlse  loading,  c'C^/c'C^  vs  X span  for  a ■ 4-1/2 
degrees  (Fig.  18). 

An  analysis  of  the  resulting  lift  curves  shows  that  for  each  wing 
an  Increase  In  results  from  the  addition  of  winglets.  This  Increase, 
shown  In  Fig.  13  for  example,  coupled  with  a slight  increase  In  C^, 
allows  the  wing  to  develop  the  required  lift  at  a lower  angle  of  attack 
and  reduce  wing  parasite  drag. 

Figure  14  shows  a typical  drag  reduction  resulting  from  winglet 
addition.  This  drag  reduction  is  accompanied  by  an  increase  in  the  nose 
down  pitching  moment,  C^,  as  shown  in  Fig.  16,  and  an  Increase  in  the 
wing  root  bending  moment  coefficient  (Fig.  17).  The  pitching  moment 
curves  show  an  Improvement  in  static  longitudinal  stability,  however, 
the  increased  nose  down  moment  must  be  balanced  by  a down  load  on  the 
horizontal  tall  which  negates  part  of  the  induced  drag  reduction.  The 
increased  root  bending  moment  is  due  to  the  outboard  shift  in  wing 
loading  as  shown  in  Fig.  18. 
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INDUCED  DRAG  COEFFICIENT  *10 


Fig.  14.  Wing  III  Drag  Polars 


LIFT  COEFFICI 


1 


Fig.  15.'  vs  Cj^  for  Wing  // 1 
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Figure  15  and  the  remaining  plots  of  vs  are  approximately 
linear,  which  suggests  that  the  Induced  drag  can  be  represented  by 

Si  - Tsh  <10> 

where  e represents  the  induced  drag  efficiency  factor.  Values  of  e are 
calculated  for  each  planform  and  the  ratio,  e'/e,  where  the  prime  indi- 
cates the  value  with  wlnglets,  is  shown  as  a function  of  aspect  ratio, 
wing  sweep  angle,  and  wing  dihedral  angle  in  Fig.  19.  The  factor  e'/e 
relates  directly  to  the  percentage  reduction  in  induced  drag 


AC, 


Di 


1 - e/e’ 


'Di 


(ID 


It  should  be  noted  that  the  above  relationship  assumes  that  the  aerody- 
namic improvement  due  to  wlnglet  addition  is  solely  accounted  for  in 
the  efficiency  factor,  e',  and  the  aspect  ratio  remains  constant. 

In  Fig.  19,  the  factor  e'/e  is  seen  to  improve  as  aspect  ratio 
increases.  This  indicates  that  the  percentage  induced  drag  reduction 
is  greatest  on  the  higher  aspect  ratio  wings.  The  percentage  drag 
reduction  increases  almost  linearly  with  increasing  wing  sweep  angle, 
but  increases  only  slightly  with  increasing  dihedral  angle. 


Longitudinal  Results 

Figure  20  illustrates  the  effect  of  wing  sweep  angle  on  induced 
drag  reduction  and  pitching  moment  increment.  For  each  sweep  angle  the 
change  in  induced  drag  reduction  is  proportional  to  the  square  of  the 
lift  coefficient  as  indicated  by  Fq  (10).  As  wing  sweep  angle  increases, 
induced  drag  reduction  increases  slightly  more  than  linearly,  indicating 
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Efficiency  Factor  vs  Aspect  Ratio,  Sweep  Angle  and  Dihedral 


that  vinglet  effectiveness  is  greatest  on  higher  swept  wings.  The 
aagnltude  of  induced  drag  reduction,  as  well  as  the  percentage.  In- 
creases with  Increasing  wing  sweep  angle.  As  wing  sweep  angle 
Increases,  the  Incremental  change  in  pitching  moment  Increases.  This 
change  is  approximately  linear  with  lift  coefficient  as  would  be 
expected  from  the  linear  moment  curves  in  Appendix  B. 

Figure  21  illustrates  the  Induced  drag  reduction  and  change  In 
pitching  moment  coefficient  as  a function  of  aspect  ratio.  Again,  the 
induced  drag  reduction  is  proportional  to  the  square  of  lift  coefficient, 
and  ACm  is  linear  with  lift  coefficient.  This  figure  shows  that  the 
magnitude  of  the  induced  drag  reduction  decreases  as  the  aspect  ratio 
increases.  Thus,  while  the  percentage  drag  reduction  is  greater  for 
higher  aspect  ratio  wings,  the  magnitude  of  the  drag  reduction  is 
greater  for  lower  aspect  ratio  wings. 

The  results  of  varying  the  dihedral  angle  (Fig.  22)  show  that 
induced  drag  reduction  improves  slightly  with  increasing  dihedral  angle. 
This  improvement  is  very  slight  as  the  angle  is  increased  from  0°  to  10* 
dihedral,  much  less  than  the  increase  from  -10  to  0*  dihedral.  This 
Figure  indicates  that  winglets  are  more  effective  when  used  on  wings 
with  significant  dihedral  and  this  is  consistent  with  Boeing's  specu- 
lation on  the  increased  winglet  effectiveness  on  the  KC-135  over  the 
C-141  (Ref  4). 

Again,  accompanying  this  increased  drag  reduction  is  a negative 
increment  in  pitching  moment. 

While  the  emphasis  should  be  placed  primarily  on  the  trends 
mentioned  above,  the  actual  magnitude  of  the  drag  reduction  due  to 
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wlnglets  is  believed  to  be  quite  good.  NASA  wind  tunnel  tests  on  the 
KC-135  (AR  - 7.035,  * 35  degrees,  X ■’  .33,  T ■ 7 degrees)  show 

that  the  reduction  in  drag  coefficient  at  Mach  .7,  ■ .4  is  AC^  » 

.0015  (Ref  5:59).  Figure  20  shows  that  for  AR  » 7,  A^y^ c ■ 35  degrees, 
and  X * 1/3  the  change  in  Induced  drag  coefficient  is  AC^  ■ .00165. 
Figure  22  suggests  that  the  effect  of  positive  dihedral  is  to  Increase 
this  figure  slightly.  It  should  be  noted  that  the  BAC  311  airfoil 
sections  used  in  this  study  is  one  of  the  sections  used  in  the  KC-135. 
Although  the  value  of  drag  reduction  calculated  in  this  study  appears 
high,  it  represents  only  the  change  in  induced  drag.  The  change  in 
total  drag  coefficient  would  be  expected  to  be  less  due  to  Increased 

viscous  drag.  The  value  calculated  in  this  study  represents  a valid 

♦ « . 

estimate  of  the  expected  induced  drag  reduction  in  this  case. 

Lateral  Results 

The  effects  of  wing  sweep,  aspect  ratio,  and  dihedral  angle  on 
the  wing  side  force  coefficient,  C^,  root  bending  moment  coefficient 
C^,  and  yawing  moment  coefficient,  C^,  are  shown  in  Figures  23-25. 

These  coefficients  were  calculated  for  the  right  half  of  the  configura- 
tion only.  The  total  values  for  symmetric  configurations  are  zero.  For 
each  wing,  the  addition  of  winglets  Increases  the  inward  side  force  on 
the  wing,  increases  the  wing  root  bending  moment,  and  decreases  the 
yawing  moment. 

In  Figure  23,  AC^  is  seen  to  increase  as  wing  sweep  angle  Increases, 
indicating  greater  inward  loads  are  produced  on  higher  swept  wings.  AC^ 
is  seen  to  decrease  with  increasing  wing  sweep  indicating  the  outboard 
shift  in  loading  is  not  as  great  on  higher  swept  wings.  AC  is  caused 
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by  the  drag  reduction  and  Increases  with  increasing'  wing  sweep.  ..dC^ 

'2 

is  also  proportional  to  C£  as  is  AC^. 

Figure  24  shows  that  ACy  increases  with  increasing  aspect  ratio. 
Although  the  inward  force  coefficient  at  the  winglet  is  greater  at  low 
aspect  ratios  (greater  lift-induced  inflow),  the  ratio  of  winglet  area 
to  wing  area  increases  as  aspect  ratio  Increases  (wing  area  decreases). 
Thus,  the  effect  of  the  winglet  is  greater  at  high  aspect  ratios.  ACR 
Increases  with  increasing  aspect  ratio  for  this  same  reason.  AC^ 
decreases  as  aspect  ratio  increases  and  AC^  decreases. 

Figure  25  shows  that  ACy  is  greatest  on  the  wing  with  anhedral 
while  ACR  and  ACN  are  greatest  for  the  wing  with  dihedral.  The  signi- 
ficant change  in  ACR  is  due  to  the  effective  elevation  of  the  winglet 
above  the  reference  axis. 


Winglet  Cant  Angle 

Appendix  C contains  comparisons  of  aerodynamic  data  for  both  wings 

with  + 3 degrees  winglet  Incidence.  For  both  wings  the  +3  degree  cant 

angle  yields  the  larger  Increase  in  C^.  Examination  of  the  drag  polars 

shows  that  for  both  wings  the  +3  degree  cant  angle  gives  the  lower 

• . 

induced  drag  over  the  greatest  range  of  lift  coefficients.  Since  the 
goal  of  winglet  addition  is  the  reduction  of  drag  at  cruise  conditions 
(C^  > .4),  it  appears  as  though  the  positive  cant  angle  would  be  the 
better  choice.  NASA  wind  tunnel  results  indicate  the  optimum  winglet 
cant  angle  is  significantly  less  than  that  predicted  by  the  theory 
used  here.  This  difference  is  due  to  viscous  effects  not  considered 
in  this  model. 
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Winslet  Penalties 

In  addition  to  the  desired  effect  of  wlnglets,  the  reduction  of 
drag,  there  are  certain  penalties  that  must  be  accounted  for.  The 
addition  of  wlnglets  adds  profile  and  skin  friction  drag  which  were 
not  accounted  for  in  this  analysis.  The  weight  of  the  wlnglets  end 
structural  modifications  requires  that  the  wing  develop  increased 
lift.  Initial  Boeing  estimates  indicate  that  the  addition  of  wlnglets 
to  the  KC-135  will  add  approximately  600  pounds  to  the  aircraft  basic 
weight  (Ref  A: 51).  There  is  also  a trim  penalty  due  to  the  nose  down 
pitching  moment  increase  and  a possible  structural  penalty  due  to  the 
Increased  wing  root  bending  moment.  These  effects  will  be  examined 
separately. 

Trim  Effects 

Not  only  do  wlnglets  reduce  the  induced  drag  on  a wing,  they 
increase  the  nose  down  pitching  moment  which  in  turn  must  be  balanced 
by  a reduction  in  lift  on  the  horizontal  tail.  To  evaluate  this  trim 
penalty  on  a typical  aircraft,  the  following  conditions  were  assumed 


h 

b 

c 

W 


A1/4c 


30.000  ft 
100  ft 
21.43  ft 

130.000  lbs 
35 


M = .7 

AR  = 7 

S - 1428.6  ft2 
CL  “ *4 

q = 227.8  lbf/ft2 


This  aircraft  is  approximately  equal  in  size,  weight,  and  characteristics 
to  the  Boeing  727.  From  Figure  B-3  it  is  found  that  the  addition  of 
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winglets  reduces  the  Induced  drag  coefficient  by  .0018.  The  static 
margin  Is  assumed  negligible  In  comparison  with  the  distance  between 
the  center  of  gravity  and  the  tall  aerodynamic  center,  By  aummlng 
moments  about  the  center  of  gravity,  the  Increased  lift  coefficient 
required  on  the  main  wing  is  found  to  be 

4cl  • « ‘Vxt  <“> 

For  this  hypothetical  aircraft,  xfc  ■ 60  ft.  From  Figure  B-29,  the 
change  In  pitching  moment  is  .0175.  Therefore, 

ACl  - .36  AC^  (13) 

and  the  required  lift  coefficient  would  now  be  - .406.  This  In  turn 
lessens  the  previous  induced  drag  reduction.  The  effect,  however.  Is 
very  slight. 

Structural  Weight  Increase 

At  the  same  flight  condition  examined  above.  Fig.  B-38  shows  that 
the  wing  root  bending  moment  coefficient  Is  increased  from  .043  to 
.045,  a 4.7Z  Increase.  This  increase  may  require  an  Increase  In  struc- 
tural weight,  particularly  in  the  outboard  regions  of  the  wing. 

As  a rough  approximation  of  the  weight  increase  required,  it  is 
assumed  that 


At 

T" 


(14) 


or  that  the  structural  skin  must  be  thickened  by  the  percentage  Increase 
in  root  bending  moment.  If  it  is  assumed  that  the  structural  box  is 
approximately  .4  c in  length,  then  the  total  area  to  be  thickened  Is  .8  S 


il 
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Assuming  the  original  material  is  1/4  inch  aluminum,  then 

ACR 

AW  - p ,8St-r 


(15) 
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In  this  example,  the  weight  Increase  Is  189  pounds. 
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V.  Conclusions  and  Recommendations 
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Wing  Parameter  Study 

In  the  preceding  chapter,  it  was  pointed  out  that  the  induced  drag 
reduction  resulting  £rom  wlnglet  addition  was  greatest  for  wings  with 
high  wing  sweep,  low  aspect  ratio,  and  positive  dihedral. 

The  percentage  reduction  in  induced  drag  was  greatest  for  high  I 
aspect  ratio,  and  positive  dihedral  wings.  r 

Accompanying  the  reduction  in  induced  drag  in  each  case  was  a 
negative  increment  in  pitching  moment  and  an  increase  in  the  wing  root 
bending  moment.  Both  of  these  moment  increments  negate  some  of  the  in- 
duced drag  reduction,  but  the  effects  of  each  were  very  slight.  Induced 
drag  reductions  in  the  range  of  15%  to  25%  appear  possible  for  the  wings 
considered  in  this  study.  A limitation  of  this  result  was  the  inability 
of  the  model  to  predict  viscous  and  aeroelastic  effects. 

This  study  has  indicated  that  significant  drag  reductions  are 
possible  on  wings  other  than  the  transport  wings  studied  to  date.  The 
curves  shown  in  the  previous  chapter  may  be  used  as  a first  approxima- 
tion to  the  induced  drag  reduction  possible  for  a variety  of  wings. 

Winglet  Cant  Angle  Study 

/ 

The  variation  of  winglet  cant  angle  in  this  study  indicates  that 
positive  cant  angles  produce  the  larger  inward  side  forces  on  the  wing- 
lets  and  offer  better  drag  reduction  over  the  greater  range  of  lift 
coefficients.  This  increased  loading,  however,  may  be  undesirable.  It 
should  be  noted  that  the  results  of  this  theoretical  analysis  do  not 
agree  with  the  KAf’A  theory  and  experimental  results  which  included  vis- 
cous (.fleets  not  considered  here. 
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Recommendations 

Although  the  NFLS  theory  and  method  of  calculating  aerodynamic 
coefficients  are  very  useful  and  have  yielded  useful  results,  the 
calculations  are  only  theoretical  predictions  and  should  be  considered 
as  such.  Future  investigations  should  couple  wind  tunnel  results  with 
theoretical  predictions.  Future  investigations  to  experimentally  find 
an  optimum  winglet  cant  angle  would  be  very  beneficial. 


Bibliography 


1.  Klineberg,  J.  M. , Kramer,  J.  J.,  and  Povinelli,  F.  P.  "Improving 
Aircraft  Energy  Efficiency",  Astronautics  and  Aeronautics. 

14:18-31  (February  1976).  ————————————— 

2.  Jacobs,  Peter  F.,  and  Flechner,  Stuart  G.  The  Effect  of  Winglets 
on  the  Static  Aerodynamic  Stability  of  a Representative  Second 
Generation  Jet  Transport  Wing"  NASA  TN  D-8267,  Washington: 

National  Aeronautics  and  Space  Administration,  July  1976. 

3.  Whitcomb,  Richard  T.  A Design  Approach  and  Selected  Wind-T|"i,|f^ 

Results  at  High  Subsonic  Speed  for  Wing-Tip  Mounted  Winglets.  • 
NASA  TN  D-8260,  Washington:  National  Aeronautics  and  Space 

Administration,  July  1976. 

4.  Ishimitsu,  K.  K.  Design  and  Analysis  of  Winglets  for  Military 

Aircraft.  Seattle:  Boeing  Commercial  Airplane  Company, 

February  1976. 

I 

5.  Goldhammer,  M.  I.  A Nonplanar  Lifting  Systems  Method.  Long  Beach, 

California:  McDonnell  Douglas  Corporation,  January  1976. 

6-  • A Lifting  Surface  Theory  for  the  Analysis  of  Nonplanar 

Lifting  Systems.  AIAA  Paper  No.  76-16,  Washington:  American 

Institute  of  Astronautics  and  Aeronautics,  January  1976. 

7.  Riley,  D.  R.  Wind  Tunnel  Investigation  and  Analysis  of  the  Effects 

of  End  Plates  on  the  Aerodynamic  Characteristics  of  an  Unswept  Wing. 
NASA  TN  2440,  Washington:  National  Aeronautics  and  Space  Administra- 

tion, August  1951. 

8.  Roskam,  Jan.  Flight  Dynamics  of  Rigid  and  Elastic  Airplanes. 

Lawrence,  Kansas:  Roskam  Aviation  and  Engineering  Corporation,  1972. 

9.  Karamchetti,  Krishnamurty.  Principles  of  Ideal-Fluid  Aerodynamics. 

New  York:  John  Wiley  and  Sons,  1966. 


y 


Appendix  A 

Nonplanar  Lifting  Systems  Theory 

A new  nonplanar  lifting  surface  theory  was  developed  by  Mr.'  M,  I. 
Goldhammer  of  Douglas  Aircraft  Company.  This  method  represents  non- 
planar lifting  systems  by  distributed  vortlclty.  It  has  been  used  to 
compute  the  Induced  drag  benefit  of  winglets  and  the  agreement  with 
NASA  experimental  data  Is  excellent  (Ref  5). 

The  Nonplanar  Lifting  Systems  (NPLS)  Method  is  a nonlinear,  non- 
planar, lifting  surface  theory;  wing  thickness  effects  are  neglected. 

This  theory  makes  no  small  perturbation  assumptions  as  do  previous 

V i 

lifting  surface  theories.  \ . 

The  basic  problem  is  to  determine  aerodynamic  forces  and  moments 
on  an  arbitrary  three  dimensional  lifting  system.  The  system  is  assumed 
to  be  placed  in  a uniform,  inviscid,  incompressible,  irrotational  flow. 
This  method  is  restricted  to  the  analyses  of  wing-like  bodies  only. 
Fuselages,  nacelles,  and  other  non-lifting  bodies  are  excluded.  . 

The  primary  assumption  in  this  theory  is  that  thickness  effects  are 
small  and  can  be  neglected.  Experience  has  shown  that  thin  wing  theories 
give  good  solutions  in  terms  of  sectional  and  total  aerodynamic  forces 
and  moments.  It  is  also  assumed  that  the  trailing  vortex  does  not  bend 
aft  of  the  trailing  edge.  Finally,  it  is  assumed  that  the  entire  flow 
field  is  subcritical. 

With  these  assumptions,  the  problem  can  be  treated  using  potential 
flow  techniques  and  the  solution  is  governed  by  Laplace's  equation 

A2*  = 0 ‘ (A-l) 


1 


- 
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This  equation  must  be  solved  subject  to  boundary  conditions  which  specify 
that  the  flow  be  tangent  at  all  points  on  the  surfaces  and  that  there  be 

no  disturbance  infinitely  far  from  the  body. 

-►  *►  ♦ 

[(U,  + u)  i+  Vj+  (W.  + w)  k)  • n - 0 

or  <A-2) 

4<«.  + ♦>*<> 


on  the  body  and. 


♦ - 0 


(A-3) 


infinitely  far  from  the  body. 

In  Eq  (A-l)  ♦ represents  the  potential  function.  In  Eq  (2), 

(u,  v,  w)  are  the  perturbation  velocity  components,  ♦ represents  the 

perturbation  velocity  potential.  U,,  W-f  and  *,  are  the  corresponding 

"*  * • • * 

free  stream  quantities  and  n is  the  unit  normal  vector  on  the  surface. 

The  solution  of  Eq  (A-l)  subject  to  the  boundary  conditions, 

Eqs  (A-2  and  A-3),  is  obtained  by  application  of  Green’s  Theorem  to 
determine  4>  (Ref  5). 

i-  - |4]  ✓ (x-O2  + (y-n)2+(z-c)2 

8n  3n 

(A-4) 

In  Eq  (A-4)  the  region  S is  the  surface  of  the  body  and  (C,n,C) 
are  dummy  (x,  y,  z)  variables  on  the  surface. 

To  integrate  the  resulting  induced  velocity  equations,  (Ref  6:3), 
three  simplifying  assumptions  are  made.  By  allowing  each  finite  element 
to  be  planar  the  integration  can  be  accomplished  in  a local  coordinate 


*(x.  y*  z)  * h 

i 
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system.  Second,  It  Is  assumed  that  y,  the  vortex  strength,  does  not 
vary  spanwise.  This  requires  a finer  finite  element  mesh  where  large 
changes  in  spanwise  loading  occur.  Third,  finite  elements  are  rectan- 
gular. An  Elementary  Vortex  Distribution  technique  assigns  a functional 
form  to  the  chordwise  vorticlty  distribution  over  each  finite  element. 
Thus,  the  vorticlty,  y,  in  the  Induced  flow  equations  can  be  replaced 
by  an  unknown  vortex  strength,  y',  times  a function  of  only  geometric 
parameters.  The  integrands  of  these  equations  now  contain  only  geometric 
parameters  since  the  vortex  strength,  y',  can  be  taken  outside  the 
integrals.  The  induced  flow  can  be  expressed  as 


q-y'(A  i + A j+A  k)+U 

’ U V J w 


(A-5) 


where  A's  are  influence  functions.  The  finite  element  solution  can  now 
be  formulated. 

Velocities  induced  by  all  finite  elements  are  summed  at  each 
boundary  condition  point.  The  vortex  strength,  y’,  is  then  adjusted 
until  the  normal  velocity  at  each  point  is  zero. 


Elementary  Vortex  Distribution 

The  Elementary  Vortex  Distribution  technique  assigns  a functional 
form  to  the  vorticlty  distribution  on  each  element.  Thin  airfoil 
theory  (Ref  9)  shows  that  the  vorticlty  distribution  on  a flat  plate 
airfoil  is  of  the  form 


Y(C)  - J C/c  (1-C/c) 


(A- 6) 


At  the  leading  edge  (£  -*■  0)  the  function  behaves  as  1/^K.  1//^,  then, 

appears  to  be  a good  choice  for  the  lending  edge  element. 
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The  actdal  form  used  In  the  present  method  Is 

Y<«  - | T'  t(C/6)1/2  - «/«>],  0 < K < « (A-7) 

The  fector  of  2/3  allows  y*  to  represent  the  mesn  value  of  over  the 
Interval.  The  leading  edge  EVD  is  shown  in  Fig.  A-l. 

Aft  of  the  leading  edge  element , y is  represented  by  a trapezoidal 
vorticlty  distribution  which  results  from  the  superposition  of  triangu- 
lar distributions. 

Y’/«l  (Cx  +«!>-«!<  ^ « 0 
Y(«)  ■ < (A-8) 

-y'/«2  u2  - «2)  o < e2  1 «2 

k 

This  distribution  is  also  illustrated  in  Fig.  A-l. 

It  should  be  noted  that  this  distribution  does  satisfy  the  Kutta 
condition  that  flow  leave  the  trailing  edge  tangential  to  the  camberline 
and  there  be  no  loading  at  the  trailing  edge. 


Calculation  of  Forces  and  Moments 

After  the  values  of  the  influence  functions  have  been  obtained,  it 
is  necessary  to  calculate  surface  pressures,  forces,  and  moments.  The 
basic  relation  used  to  compute  the  pressure  jump  across  the  lifting  sur- 
face is  the  Kutta-Joukowski  Law.  It  states  that  the  force  per  unit  length 
on  a vortex  in  a flow  field  is  proportional  to  the  fluid  density  times  the 
vortex  strength  times  the  velocity.  Aerodynamic  forces  and  moments  are 
computed  by  Integrating  the  pressure  jump  coefficients  over  the  lifting 
surfaces.  Although  the  flow  field  was  initially  assumed  to  be  incompres- 
sible, the  progrsm  allows  for  corrections  based  on  the  Prandtl  - Glauert 
transformation  to  Include  compressibility  effects. 
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Appendix  B 

Aerodynamic  Data  for  Wing  Parameter  Study 

This  Appendix  contains  plots  of  the  aerodynamic  coefficients 
calculated  for  the  wing  parameter  study  using  the  NPLS  Computer  Pro- 
gram. 

The  following  figures  show  comparisons  of  aerodynamic  data  for 

I 

the  wing  alone  with  that  for  the  wing-winglet  combination.  All  wings 
were  constructed  using  a Boeing  BAC  311  airfoil  and  had  a taper  ratio 
of  1/3.  The  winglets  used  were  desribed  on  page  and  were  mounted 
at  0 degrees  cant  angle  and  70  degrees  dihedral  angle. 

The  aerodynamic  data  presented  was  calculated  for  the  following 
wings: 


Wing  # 

A1/4c  = 

AR  = 

r= 

2 

25 

7 

0 

3 

35 

7 

0 

4 

45 

7 

0 

5 

15 

6 

0 

6 

15 

5 

0 

7 

15 

4 

0 

8 

45 

4 

0 

9 

15 

7 

+ 10 

10 

15 

7 

-10 
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Fig.  B-l.  Wing  |2  Lift  Curve*  Fig.  B-2.  Wing  IS  Lift  Curve* 
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